Constraining chameleon field theories using the GammeV afterglow experiments 
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The GammeV experiment has constrained the couplings of chameleon scalar fields to matter and 
photons. Here we present a detailed calculation of the chameleon afterglow rate underlying these 
constraints. The dependence of GammeV constraints on various assumptions in the calculation 
is studied. We discuss GammeV-CHASE, a second-generation GammeV experiment, which will 
improve upon GammeV in several major ways. Using our calculation of the chameleon afterglow 
rate, we forecast model-independent constraints achievable by GammeV-CHASE. We then apply 
these constraints to a variety of chameleon models, including quartic chameleons and chameleon 
dark energy models. The new experiment will be able to probe a large region of parameter space 
that is beyond the reach of current tests, such as fifth force searches, constraints on the dimming of 
distant astrophysical objects, and bounds on the variation of the fine structure constant. 



I. INTRODUCTION 

Increasingly strong evidence has emerged over the past 
decade that the expansion of the universe is accelerating, 
a phenomenon which can be explained by a scalar field 
"dark energy" with negative pressure [E S H, IE B 13 • 
Couplings between such a scalar and Standard Model 
particles could lead to fifth forces observable at labo- 
ratory or solar system scales, that must be hidden in 
order for the dark energy to satisfy constraints on the 
non-observation of such forces. There are three known 
ways to hide dark energy-mediated fifth forces: weak cou- 
plings between dark energy and matter d, [l(| ; effec- 
tively weak couplings locally, as in the Dvali-Gabadadze- 
Porrati (DGP) brane world model and its generaliza- 
tions [Tl]> [lH, [l3| ; and an effectively large mass locally, as 
in chameleon theories [Til [TH [l6| . 

The chameleon mechanism makes the effective mass 
of a scalar field grow substantially as the density of 
the surrounding matter is increased. Chameleon theo- 
ries are a particularly well-studied class of dark energy 

models Ha m m, m m, m m m m m. si nce 

fifth forces at least as strong as gravity are already ruled 
out from solar system scales to submillimeter scales, the 
chameleon mechanism must raise the effective field mass 
above about (0.1 mm) -1 ~ 1CP 3 eV locally. In addi- 
tion to chameleon dark energies, the "scalaron" field, 
resulting from the transformation to Einstein frame of 
an f(R) gravity theory, is able to satisfy local and stel- 
lar tests of g ravity only through the chameleon mecha- 
nism dJSiii,^. 

It was shown in [3l|, that chameleons which couple 
strongly to matter have correspondingly large chameleon 
effects, allowing them to evade laboratory fifth force con- 
straints such as [33] on matter couplings (3 m ~ 1. If these 
chameleons also couple strongly to photons, then they 
are ideally suited to probes of the afterglow phenomenon, 
the first of which was the GammeV experiment at Fermi- 
lab [l7|. Consider a closed, evacuated cylindrical cham- 



ber, with glass windows at the ends and a magnetic field 
in the interior. Photons streamed through the windows 
will occasionally oscillate into chameleon particles in the 
background magnetic field. If the mass of one of these 
chameleons in the walls of the chamber is greater than its 
total energy inside the chamber, then it will reflect from 
the wall; such a chameleon will be trapped in the cham- 
ber. After the photon source has been turned off, any 
remaining chameleons will oscillate back into photons in 
the magnetic field, producing an observable "afterglow" 
of photons. 

Here, we study the oscillation between photons and 
chameleon particles in the cylindrical vacuum chambers 
used in afterglow experiments such as GammeV. We com- 
pute the afterglow and decay rates per chameleon particle 
as a function of the dimensionless chameleon-photon cou- 
pling constant /3 7 . Our calculation allows for absorptivity 
in the chamber walls, chameleon-photon phase shifts due 
to a nonzero chameleon mass m e s (chamber), and addi- 
tional phase shifts fref due to reflections from chamber 
walls, as found in [18j . For a particle in a superposition 
of chameleon and photon states, the glass windows at 
the ends of the chamber act as quantum measurement 
devices; photons are transmitted, while chameleons are 
reflected. The chameleon decay and afterglow rates can 
be computed by keeping track of the decline in chameleon 
amplitude between successive quantum measurements by 
the windows. 

Armed with these rates, we proceed to compute the ex- 
pected afterglow flux as a function of time in GammeV, 
as well as the resulting constraints on the chameleon- 
photon coupling as a function of the chameleon mass, a 
result that was discussed in ref. [i~7j] . We then forecast 
the constraints that will be attained by an upcoming 
experiment, GammeV-CHASE, expected to take data 
in the winter of 2009-2010. GammeV-CHASE was de- 
signed to improve in several ways upon the first GammeV 
chameleon experiment, hereafter referred to as GammeV. 
We show that GammeV-CHASE will improve constraints 
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on large photon couplings /3 7 by several orders of mag- 
nitude, bridging the gap between GammeV constraints 
and the collider constraints of ref. [l9j . Furthermore, it 
will reach chameleon masses several times greater than 
those excluded by GammeV, allowing us to probe masses 
at the dark energy scale p^{ 4 = 2.4 x 1CV 3 eV. 

Finally, we use our afterglow computation to fore- 
cast GammeV-CHASE constraints on specific chameleon 
models. We begin with power law potentials, de- 
voting most of our attention to </> 4 potentials, which 
are well-understood and ubiquitous in particle physics. 
GammeV-CHASE constraints on cj) 4 chameleons will 
complement the fifth force experiments of Eot-Wash [33[ 
by probing strong matter couplings. Next, we study 
two different models of chameleon dark energy, with in- 
verse power law potentials and exponential potentials. 
GammeV-CHASE will probe large areas of the f3 m , /3 7 
parameter space of each of these dark energy models. 
These parameter regions are largely unexplored by the 
current data, including constraints on the dark energy 
equation of state, variations in the electromagnetic fine 
structure constant, and the dimming of distant objects 
through photon-scalar oscillation in the Galactic mag- 
netic field. Thus, the GammeV-CHASE laboratory 
search for afterglow will complement current cosmolog- 
ical probes of dark energy, vastly extending our con- 
straints on dark energy couplings to matter and photons. 

The paper is organized as follows. After a brief 
discussion of chameleon physics in Sec. [Til we study 
chameleon-photon conversion in afterglow experiments 
in Sec. IIII1 and compute the decay and afterglow rates 
for GammeV and GammeV-CHASE in Sec. UVl We 
constrain chameleons using GammeV and forecast con- 
straints for GammeV-CHASE in Sec.fV] before conclud- 
ing in Sec. |VU 



II. CHAMELEON FIELD THEORIES 



density p m , 



A. Equations of motion 



We consider chameleons with action 
S = j d i x^(^Ml l R-^d li 4>d^4>-V{4>) 



2/3 m <p/M P 



9n»,4>m)) (!) 



where /3 7 and (3 m are, respectively, the dimensionless 
chameleon couplings to photons and matter. We note 
that elsewhere in the literature, these couplings are ex- 
pessed in dimensional form as <? 7 = /3 7 /Mpi and g m = 
Pm/Mpi: for a typical value (i 1 = 10 12 probed by Gam- 
meV [l7j], we have g 1 = 4.1 x 1(T 7 GeV -1 . £ m is the 
matter Lagrangian; the matter coupling (3 m is assumed 
to be universal to all species of matter. Varying with 
respect to 6 an( i the electromagnetic field, we find the 
equations of motion in the presence of a constant matter 



06 = - 



dV cB 



(2) 



Vcff {x, 6) = V{6) + e%% m + ie^pt F^F^ (3) 
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with the other two of Maxwell's equations unchanged. 
That is, the field couples to the matter density p m 
and the electromagnetic field Lagrangian density p 1 = 
(F^F^/A = {\B\ 2 - \E\ 2 )/2. Perturbing about a con- 
stant background scalar (f> and magnetic field B, we see 
that oscillations from chameleon to photon are described 
by 
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V x (B x V)S6 = 
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to first order in the perturbations [3J 



B. Chameleon effects 

The field value (f> m m at the minimum of the effec- 
tive potential, and, hence, the effective mass m 2 s = 
V off . 00(<^min); can vary with the background matter den- 
sity and electromagnetic fields. If m G fj increases signifi- 
cantly with p m or p 7 , then the field is called a chameleon 
field. By acquiring a large mass at laboratory and so- 
lar system densities, the chameleon can "hide" from 
fifth force constraints on ordinary matter-coupled scalars. 
As a simple example, consider a power law chameleon, 
V{4>) — gS^ ' i in a background with negligible E and B. 
Assuming that /3 m m i n /Mpi <C 1, which is required by 
current constraints, the mass scales as m c ff oc with 
rj = [M - 2)/{2M - 2). If M = 2, that is, the potential 
is just a mass term, then i] = 0; m c g does not scale with 
density, and the field is not a chameleon. For Af = 4, 
on the other hand, 77 = 1/3, and references [H], H2] show 
that, for sufficiently large /3 m , chameleon effects allow this 
theory to evade fifth force constraints from torsion pen- 
dulum experiments such as Eot-Wash (33|. Chameleon 
fields strongly coupled to matter acquire large masses 
inside the layer of foil used to keep the source and test 
masses electrostatically isolated. If the chameleon Comp- 
ton wavelength m~a is much less than the thickness of 
the foil, then chameleon effects in the foil will screen the 
fifth force of the source mass on the test mass, rendering 
the chameleon field undetectable. 



III. CHAMELEON-PHOTON CONVERSION IN 
AFTERGLOW EXPERIMENTS 

A. Photon production 

Let the chameleon and photon amplitudes be = 5cf> 
and ^ 7 = SB/k, respectively. We consider the oscillation 
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into photons of chameleons that are already trapped in- 
side the GammeV chamber for the case of a small mixing 
angle. Thus we assume that |\& 7 | = initially, and that 



|\1/ 7 | <C |^0| ~ 1 throughout the calculation. In the 5 



presence of a constant magnetic field in the x direction, 
B = Bx, (O becomes 



2L 

dt 2 



- k' 



*7 = 



hh^-k x (x x k)^4,. 



(6) 



In the relativistic limit, {—d 2 /dt 2 — \k\ 2 ) — (id/dt 
k)(id/dt + k) sa 2k(id/dt ~ k), this has the solution 



# 7 (t) = 



-ikt t J l 



(3 7 B 



a(k) exp 



—ikt — 



2Mpi jq 

i 2kf3-yB 



'off 
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sin 



2k 



dt 7 (7) 



4A- 
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where m c ff is the effective mass of the chameleon, 
k = kk = fc(sin(0) cos(f)x + sin(0) sin(^)y + cos(#)z) 
is the particle momentum, and a(k) = k x (x x 
fc) = (1 — sin 2 (#) cos 2 (y))i — sin 2 (#) sin(<^) cos(<p)j/ — 
sin(0) cos(9) cos((p)z. Using \a\ 2 — sin 2 (<,5) + 
cos 2 (6) cos 2 (ip), and setting 9 = 0, we recover the fa- 
miliar formula for the oscillation probability for a path 
perpendicular to B, 



C 2 sin 2 



"■off' 

4fc 



(9) 



where we have defined C = 2k(3 1 B / (Mpxm 2 ^) . 

Next, assume that the chameleon particle is inci- 
dent upon a wall, with normal vector h. We consider 
chameleon models for which the effective mass inside 
the wall is much greater than the chameleon energy 
u> = \Jk 2 + m 2 ff inside the chamber. Thus the chameleon 
cannot penetrate the wall, and tunneling is negligible; the 
particle must bounce. The bounce has three effects: 



1. the direction changes, k — > k — 2(k ■ h)n, with a 
corresponding change in a; 

2. the photon component has a probability / a b s = 1 — 
/ re f of being absorbed in the walls; 

3. the photon is phase shifted by an angle £ rc f relative 
to the chameleon UM. 



These last two imply that ^ 
where A rcf = /^ 2 e ,&ef , 



A-cf ^7 due to a bounce, 



B. Afterglow experiments 

A simple afterglow experiment will trap chameleon 
particles in a cylindrical chamber, of radius R, with glass 
windows at the entrance and exit. A magnetic field re- 
gion of length L inside the cylinder will be offset from the 
entrance by a length and from the exit by a length £2, 




FIG. 1: Particle in an afterglow experiment. The magnetic 
field region is shaded. Inside this region, where chameleon- 
photon oscillation takes place, the particle trajectory is drawn 
as a dashed path. Next to each dashed segment is shown the 
contribution of that segment to the total photon amplitude 
for the case where / re f = 1, £ re f = 0, and m ef r <g 4nk/L. 



for a total chamber length of itot = £1 + £2 + L. Photons 
stream into the chamber through the entrance window, 
and out through the exit, oscillating into chameleon par- 
ticles with a probability given by ©. For a particle in 
a superposition of chameleon and photon states, a glass 
window serves as a quantum measurement device; a pho- 
ton will pass through, while a chameleon will bounce. 
Thus we can assume that any particle bouncing from ei- 
ther window is in a pure chameleon state. After the pho- 
ton source is turned off, a photomultiplier tube (PMT) 
in a dark box outside the exit window will begin to look 
for afterglow photons produced by chameleons trapped 
in the chamber. 

Let the origin of the coordinate system be at the center 
of the entrance window, with the cylinder extending in 
the positive z direction, and the magnetic field in the x 
direction. Assume that the particle begins at the origin; 
we will test the effects of this assumption later. Let 9 
and tp be defined in the usual way; 9 is the angle between 
k and the z axis, and ip is the angle between k — (k ■ 
z)z and the x axis, where k is the initial direction of 
the particle. For small 9, the particle will not bounce 
inside the field region, and ([9]) suffices to compute the 
probability of photon production. 

At larger 9, the particle will bounce inside the magnetic 
field region. This case is sufficiently complicated that we 
begin with a simple example in which the chamber walls 
are perfectly reflective, the phase shift associated with 
wall reflection is zero, and the chameleon mass is small 
enough that phase differences between the chameleon and 
photon as they propagate through the chamber may be 
neglected. Figure [T] shows a sample particle trajectory. 
The particle begins in a pure chameleon state at the en- 
trance window, shown at the left end of the chamber, 
and heads to the right, towards the magnetic field. Let 
z\, Z2, ■ ■ ■ , Zn be the sequence of z values at which the 
particle bounces from a wall. For a particle beginning 
at the origin, z n = (n — l/2)Az, with Az = 2Rcot{9). 
Let the leftmost and rightmost bounces inside the mag- 
netic field be riL and riR, respectively; the total number 
of bounces in the B field region is then ng = tir — «l + 1 • 

The sample trajectory shown in Fig. [1] has six bounces, 
four of which are in the B field region, and is made up of 
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seven segments, five of which have some overlap with the 
B field region. In the figure, a dashed line denotes that 
portion of a trajectory which is in the B field region, and, 
hence, contributes to the photon amplitude. The contri- 
bution of each segment is given by ([8]). In our approxi- 
mation that m e ff is small, the chameleon and photon are 
in phase; "I^,, ^ 



7 oc exp(— ikt). Thus we can neglect the 



to a bounce from the walls. Even a polished metal sur- 
face will not be perfectly reflective; the fraction of in- 
cident photons reflected back into the chamber will be 
around / rc f ~ 0.9. Furthermore, since chameleons and 
photons bounce at slightly different distances from the 
wall, the bounce can introduce a nonzero phase shift 
Cref [HI- We model absorption in the walls by mul- 



each bounce. Thus the photon amplitude just before the 
HLth bounce is \& T (2~ ) = CsLa nL , and the amplitude 
just after that bounce is ^ 7 (z+ L ) = A IQ fCsi,a nL . Be- 
fore and after the next bounce, we have ^j(z, 



phase factor in front of \& 7 , and write the contribution tiplying the photon amplitude by A 
of the nth segment as C sm(m^ s t / (4k))a n , where t is the 
time spent inside the B field region, o„ = a(k n ), and k n 
is the direction of the particle before the bounce at z n . 
Furthermore, Csin(m^ ff t/(4fe)) » (3 7 Bt/(2M P1 ) at low 
m ff , so the contribution of each segment is proportional 
to its length inside the magnetic field region. 

The leftmost segment in the B field region, the seg- 
ment between bounces til — 1 and til (in Fig.[TJ bounces 
1 and 2) is only partially inside the magnetic field, so t-^ — 
(z„ L — ^i)sec(0) for that segment. Similarly, the right- 
most segment in the magnetic field, between bounces 
and riR + 1 (here, 5 and 6) , has 4r = {l\+L — z, iR ) sec(8) . 
Meanwhile, all of the segments in the middle, < n < 
tir, spend an equal amount of time in the magnetic field, 
tu — Azsec(6). Defining 



rcf 



at 



and * 7 (z+ L+1 ) = A^ ci Cs L a nij + 
A re fC 'sMflnL+ij respectively, where we have assumed that 
til, + 1 < fiR. Summing over all of the bounces, we find 
the photon amplitude at the exit of the B region, 



yjv(low mass) 
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t B ) = A£ C SL a„ L + Cs M J2 43 +1-n 3n 

n 

+Cs R a nR+ i, (13) 
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sr 



4fccos(6>) 



'cir 



Az 



4fccos((9) 
m 2 eS (ti 4 



'•cir 



R 



2ksm{9) 



4fccos((9) 



(10) 

(11) 
(12) 



we see that each of the the leftmost, middle, and 
rightmost segments contributes CsLa„ L , Csmo„, and 
CsrOmr+i, respectively, where < n < tir. The 
photon amplitude at the exit z = £g =£i+Lof the 
B field region is simply the sum of these contributions, 
* 7 (^b) = Cs L a nij +J2n Cs M a n +Cs R d nR+ i, up to a com- 
plex phase factor which disappears when \l/ 7 is squared. 

At this point, our choice of initial conditions simpli- 
fies the problem. Recall that we have assumed the initial 
particle position to be the center of the entrance win- 
dow. The trajectory of such a particle will remain in 
the plane spanned by the z axis and the initial direc- 
tion k, even after the particle bounces from the chamber 
walls. This is because the particle momentum has only 
a radial component oc xx + yy and a component propor- 
tional to z. Each bounce reverses the sign of the radial 
component while leaving the z component unchanged. 
Furthermore, each bounce reverses the sign of the i com- 
ponent of a while leaving the x and y components un- 
changed; for all odd n, a n — oi, and for all even n, 
(in — &2, with CL2 — <2i — 2(<zi • z)z. The summation 
over segments inside the B field is considerably Simpli- 
fied, J2n Cs Ma n = Cs M En( ffl t^ + a i-.vV + {-^T °-l.zZ) ■ 

Next, we consider two complications to the problem: 
the possibility of photon absorption in the chamber walls, 
and the phase shift between chameleons and photons due 



where n B — — + 1 is the total number of bounces 
in the B region. 

The final effect which we have neglected thus far is the 
phase shift associated with the chameleon-photon mass 
difference. Eq. {8]| implies that contributions to \I/ 7 from 
different segments of the path, which occur at different 
times, will have relative phases. As before, the leftmost 
and rightmost segments in the B region contribute dif- 
ferent amounts to the phase shift, since they are only 
partially inside the B region, while the middle segments 
all contribute the same phase shift. Defining 
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(14) 
(15) 
(16) 

(17) 



we have our final expression for the photon amplitude at 
the exit of the B region, 



* 7 (^b) = -iCexp 



ik£ B 
cos(<9) 



M 



+ A riB 3 s M a„ L+J ■ +e^ R s R a nR+1 



A nB e^s h a nii 
, (18) 



where we have simplified our notation by combining A rc f 
and exp(i^M) into A. 

As above, our choice of initial conditions allows us to 
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compute explicitly the sum in ([TBI) . 



"B-i A- A riB 

Y AnB ~ J ^+j = (ai, x x + a ltV y) 
3=1 

+a 1 ^(-l)- ( ' A) 1 ' + ( / ) " B .(19) 

After the particle leaves the magnetic field region, the 
photon and chameleon states decouple, and no further 
oscillation occurs. Thus, the probablity that the particle 
will be a photon, when a quantum measurement is made 

by the exit window, is |* 7 (M| 2 /ref ™ R , where N ~ n R 
is the number of times that the particle bounces after 
leaving the field region. 

Since we also want to compute the decay rate of the 
chameleon, we must keep track of the photons lost to ab- 
sorption in the chamber walls, inside the magnetic field 
region. At each bounce, \5 7 — * A re f^ T , the absorption 
probability is incremented by (l-/ rcf )|* 7 | 2 = Jabs | * 7 1 2 ; 
since the tunneling rate is negligible, the sum of the re- 
flection and absorption probabilities must be equal. We 
find the total absorption probability T^bs by summing 
Jabs | | 2 over all of the bounces in the field region. 



abs 



/absC 



3=0 



/rcfK, 



+4 Y 

3=1 

T lB— 1 j 

+SLSM Y, Ys ' 
3=1 1=1 



x (e 4 « L A> \A*y- 1 + e-^iAyAi- 1 ) 



(20) 



Once again, all of the summations can be carried out 
explicitly for particles originating at the center of the 
entrance window: 



Y fL\ 
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1 — /ref 
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3 = 1 



YA^d 
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(1-A)(1-A*) 



riB - 1 
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J rcf J lc f 
1 - /ref 
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1 + A V 1 + A* 
complex conjugate. 
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A chameleon particle with initial direction specified by 
9 and (p traverses the chamber in a time £tot sec(0). It has 
a probability V a b s + 1 (&B ) | 2 of producing a photon, and 
a probability \^^(£b) \ 2 f^f nn of producing a photon that 
escapes through the exit window. Thus the contribution 
of this 9 and ip to the decay rate is the photon production 
probability per unit time, ("P a bs + 1 * 7 {(-b ) | 2 ) / (Act sec(0) ) , 
and similarly for the afterglow rate. We find the total 
decay and afterglow rates per particle by integrating over 
angles. 
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2ti" Jo Jo 
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cos(e) 



(24) 



,N(0)-n R (e) 
/ref 



cos(0) 

Itot 
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The extra factor of two in the decay rate accounts for 
chameleons which begin at the exit window and travel 
toward the entrance window. Note that all of the depen- 
dence on ip is due to dot products of the a vectors. Since 
these differ only in the sign of the z component, they can 
all be rewritten in terms of squares of the components 
of di, as in (TiUl |2"T1 12"21 12"3"| . and the <p> integral can be 
computed explicitly using 



-!- / \t x d^ = cos 2 (9) + lsm 4 (9) 
2tt Jq a 

± J\l y dp> = \sin\9) 

±j \ 2 iZ ^ = isin 2 (0)cos 2 (0). 



(26) 
(27) 
(28) 



We need only integrate numerically over one variable, 9, 
in order to compute the decay and afterglow rates per 
particle. 

We note that the magnetic field strength B and 
the photon coupling /3 7 only appear in rdcc,7 and T a ft 
through factors of C 2 in |^ 7 | 2 and T^bs- This C 2 can be 
brought outside the integrals in (f2"4")) and ([231) . Thus the 
decay and afterglow rates scale as idcc,7,r a ft <x B 2 f3 2 , 
for any m e g , £ rc f , and chamber geometry. We need only 
compute the decay and afterglow rates for one value of 
each of B and /3 7 . 
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FIG. 2: Diagram (side view) of the apparatus used by Gam- 
meV and GammeV-CHASE to search for a chameleon after- 
glow. Diagram is not to scale. Actual values of these lengths 
for each experiment are given in Table [T] 



Next, we compute the expected flux of afterglow pho- 
tons through the exit window for this simple afterglow 
experiment. We assume that chameleon-photon oscilla- 
tion is the dominant contributor to the total chameleon 
decay rate, and that all other decays are negligible. The 
experiment proceeds in two stages: production and af- 
terglow. During the production stage, photons stream 
through the chamber at a rate of F 7 . Each photon 
has a probability V^-y = C 2 sin 2 (m^ ff L/(4fc)) of pro- 
ducing a chameleon particle, for a total production rate 
of F 7 7 , 0^ 7 . Meanwhile, if N$ chameleons are present 
in the chamber, then the total decay rate is N^T^^. 
Thus the number of chameleons is given by dN^(t)/dt 



Assuming that chameleon pro- 



duction occurs over a time interval 
have 



-r m . < t < 0, 



N 



(prod) 



(*) 



FJP, 



7 / 0<— s-7 
■ dec, 7 



l 



y (t+T pr ) 



(29) 



in the production stage. In the afterglow stage, t > 0, 
the photon rate is reduced to zero, and the number of 
chameleons is 



fVV^7 a e -r dcCiT r pr -j e -r dcc , 7 t 

Tdec.7 



(30) 



The expected signal during the afterglow stage of the 
experiment is found by multiplying by r a f t , 



F ait (t) = F -> P +^ Taa (1 

dec, 7 



(31) 



IV. GAMMEV AND GAMMEV-CHASE 
EXPERIMENTS 

A. GammeV geometry 

The actual apparatus used by GammeV to search for 
chameleons differs in a few ways from the simple experi- 
ment described above [l7j]: 

1. the GammeV chamber has a second section, with 
a smaller radius, extending from z = i\ + £ 2 + L to 
z = l\ + £2 + L + £ 3 , with the exit window at the 
end of the second section; 



quantity 


f ■' \ 7 

Gamme V 


GammeV-UHAbrj 


l\ 


2.36 m 


2.0 m 


r 

L 


6.0 m 


6.0 m 


£2 


1.16 m 


1.16 m 


n 


2.51 m 


0.30 m 





z.uo m 


n of! m 


n 

«5 


0.10 m 


0.10 m 


Ri 


2.38 cm 


3.175 cm 


R2 


1.75 cm 


3.175 cm 


^"lcns 


z.o 1 cm 


z.o 1 cm 


fPMT 


2.5 mm 


2.5 mm 


to 


1006 sec 


1 sec 


At 


OU1U ScL 


i nnn q^<^ 

^ luuu oec 


fiel 


0.53 


0.53 


k W U! 


2.33 eV 


2.33 eV 


B 


5.0 Tesla 


5.0 Tesla 


Vpump 


0.026 m 3 


not available 



TABLE I: Properties of the GammeV and GammeV-CHASE 
experiments, including the dimensions shown in Fig. [2] 



2. outside of the exit window is a lens that focuses the 
afterglow emerging from the chamber; 

3. afterglow photons must enter the aperture of the 
PMT, with radius ^pmTj m order to be detected; 

4. the PMT is uncovered a time to after the laser has 
been turned off, and data are collected for a total 
time At; 

5. the vacuum inside the chamber is maintained by 
a turbomolecular pump connected to a roughing 
pump, which increases the total volume accessible 
to chameleons by V^, U mp; 

6. chameleons light enough to enter the roughing 
pump will be removed from the chamber, so the 
range of masses which can be constrained is lim- 
ited. 

Figure[2]is a diagram of the GammeV apparatus, showing 
the lengths and radii defined above; the numerical values 
of these quantities are listed in Table|U We define £tot,i = 
£1 + £2 + L to be the length of the first section, and 
^tot,2 = 4ot,i + ^3 to be the total length of the chamber. 

Limitations imposed on GammeV by the pumping sys- 
tem were discussed in [l7j . A chameleon particle will be 
pumped out of the vacuum chamber if its effective mass 
77i ff (rough) at the intake of the roughing pump, where 
the air pressure is P r ough = 1.9 x 10~ 3 torr, is less than 
its energy u> sa k inside the chamber. This pumping oc- 
curs on time scales much smaller than to, so a chameleon 
model whose particles can be pumped out of the chamber 
is inaccessible to GammeV. In particular, we shall see in 
Sec. [V] that GammeV can only constrain a small subset 
of chameleon dark energy models. For the remainder of 
this section we assume that m e s (rough) > k. We proceed 
to calculate the decay and afterglow rates as functions of 
m e ff (chamber), /3 7 , and £ ro f, noting that any resulting 
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FIG. 3: Face, 7 (thick lines) and r a f t (thin lines) vs. m g, for 
several values of £r e f , assuming a magnetic field of B = 5 Tesla 
and a photon coupling of /3 7 = 10 12 . 
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FIG. 4: Fdcc,7 (thick lines) and F a f t (thin lines) vs. £ rc f, 
for several values of m e ff, assuming a magnetic field of B = 
5 Tesla and a photon coupling of /3 7 = 10 12 . 



constraints on chameleon models will only be applicable 
when the above mass condition is met. 

When calculating the afterglow rate, we must check 
whether the particle enters the second section, reaches 
the lens, and enters the aperture of the PMT. It is ap- 
parent that a particle beginning at the center of the en- 
trance window, our chosen initial condition, will have a 
higher probability of being near the central axis of the 
chamber, and therefore a higher probability of reach- 
ing the second section, the lens, and the PMT. Com- 
puting these probabilities for the GammeV geometry 
is straightforward; a particle beginning at an arbitrary 
position on the entrance window has a probability of 
Pavg = 2.04 x 1CT 5 of reaching the PMT, while a particle 
beginning at the center of the window has a probability 
of Pctr = 3.35 x 10~ 5 . Thus we must normalize our after- 
glow results by / ge0 m = Pavg/Pctr, which takes the value 
0.610 in this case. Henceforth, we assume that r a f t has 
been normalized appropriately. 

The final results of the computation detailed above are 
the afterglow rate r a f t per chameleon particle, and the 
rate rdec,7 of the decay of a chameleon particle to a pho- 
ton, both of which are functions of the background mag- 
netic field B = Bx, the chameleon mass m e ff (chamber) , 
the chameleon-photon coupling /3 7 , and the phase shift 
£ rc f at each wall reflection. As noted earlier, rd ec ,7 
and r a ft both scale as the square of the magnetic field 
strength and the photon coupling. The variations of these 
rates with m c g (chamber) and £ ro f are more complicated. 

Figure [3] shows Fd cc ,7 and r a f t as functions of mass 
for several values of £ re f. Consider first the thin, solid 
line, which shows r a f t for £ rc f = 0. The afterglow rate 
is dominated by chameleons on trajectories with small 9, 
which bounce no more than a few times from the cham- 
ber walls, since these are the most likely to produce pho- 
tons which reach the detector. For a trajectory with no 
bounces, © implies that the afterglow rate is propor- 
tional to sin 2 (m^ ff L/(4fc)). Thus we expect r a f t to be 



small when m cS = y/4njk/L = 9.81 x 10" 4 j 1/2 eV for 
any positive integer j; at these masses, destructive in- 
terference in the chameleon-photon oscillation suppresses 
photon production. As expected, Fig. [3J shows the first 
local minimum of the afterglow rate at m e fi ~ 0.001 eV 
and the fourth around 0.002 eV. Furthermore, these local 
minima are very sharp, since paths with slightly larger 9 
differ only slightly in phase at the end of the magnetic 
field region. In contrast, the other two afterglow rate 
plots have much shallower minima at different masses. 
This is because a nonzero £ re f implies that paths with 
different numbers of bounces can exit the magnetic field 
region with very different phases. Meanwhile, the decay 
rates, plotted as thick lines in Fig. [3J have fewer features, 
since they average paths over a much larger range of an- 
gles. 

Tdec,7 and r a f t are shown as functions of £ re f, for a few 
fixed masses, in Fig. [5] In the low-mass limit, the only 
phase shift between the chameleon and the photon is due 
to reflection from the walls. A total phase difference of 
zero implies maximal constructive interference between 
the chameleon and photon wavefunctions, and, hence, 
maximal chameleon-photon oscillation. Thus the decay 
and afterglow rates will peak at £ re f = 0. This is consis- 
tent with the plots of rd ec ,7 and r a f t for m c s = 10~ 4 eV, 
shown as thick and thin solid lines, respectively. At larger 
masses, this maximum in the decay rate shifts to nega- 
tive £ rc f, in order to compensate for the positive phases 
£l, £mj and £r caused by a nonzero chameleon-photon 
mass difference. Meanwhile, at m c g — 0.001 eV, £ rc f ~ 
corresponds to maximal destructive interference, as dis- 
cussed above. Thus a nonzero £ re f increases r a f t by mak- 
ing the interference in chameleon-photon oscillation more 
constructive, leading to a minimum around £ re f = 0, as 
shown in Fig. 2] 

Next, we will show that the effects of the initial con- 
ditions on the dynamics of chameleon-photon oscillation 
are negligible; the normalization / gcom is the only correc- 
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to [m] 


1 riar* ru H 7 

- 1 - acc,7 L J 


r a f t [Hzl 





4.905 x 10" 5 


6.643 x 10" a 


1 


4.899 x 10~ 5 


6.865 x 10" 9 


10 


4.903 x 10" 5 


6.845 x 10" 9 


100 


4.897 x 10 -5 


6.876 x 10" 9 



TABLE II: Decay and afterglow rate vs. Iq, These have been 
computed for the GammeV geometry, assuming a reflectivity 
°f /ref = 0.53, a magnetic field of B — 5 Tesla, a chameleon 
mass in the chamber of m e g = 10~ 4 eV, and a chameleon- 
photon coupling of (3-y — 10 12 . 




FIG. 5: Projection into the xy plane of a 2-point path (solid 
green line), a 3-point path (dashed red line) and a 4-point 
path (dash-dotted blue line). 



tion that we need to make to the afterglow rate. Since 
z = and \k\ are fixed, the full set of possible initial 
conditions on the entrance window is described by four 
parameters: x, y, 9, and ip. Until now, we have restricted 
ourselves to a two-parameter subset, x = y = 0. Our 
subset consists of paths that remain in the same plane as 
they bounce around, a plane which contains the z axis. 
This condition will remain true if we allow nonzero initial 
x and y, subject to the constraint y/x = tan(y). Thus 
we can study a third parameter out of the four. This can 
be done through a simple modification of our previous 
calculation. Assume that the chamber were extended by 
a length £q in the negative z direction, with the entrance 
window now at z = —£q. A particle beginning at the cen- 
ter of this new entrance window would, when it reached 
2 = 0, satisfy y/x = tan(yj), with x and y not necessarily 
zero. If we compute V a hs and |^ 7 (^s)| using these new 
paths, while using the chamber length without £q to find 
the time taken by each path, then we obtain decay and 
afterglow rates that can be compared with our previous 
calculations. Table HTl shows that varying £q, from zero 
to much larger than the length of the entire chamber, 
leaves the decay rate essentially unchanged, and changes 
the afterglow rate by no more than a few percent. 



Each path in the three-parameter subset considered 
above was confined to a single plane. Other paths do not 
remain in a plane, but travel through the chamber in a 
corkscrew fashion. It is useful to consider the projection 
of these paths onto the xy plane, as shown in Fig. [5l Let 
us also assume for this discussion that 9 is not small, 
so that each path bounces multiple times. The paths 
confined to a plane would, in this projection, appear 
to bounce back and forth between the same two points, 
passing through the center between any two bounces, as 
shown by the solid line in Fig. O We call such paths 2- 
point paths. If the condition y/x = tan(y>) were violated 
by a small amount, then the points between which the 
path bounced would appear to move slightly as the parti- 
cle progressed down the chamber; the solid line in Fig. [5] 
would appear to precess. For larger deviations, the xy 
projection of the path would close on itself, forming an 
equilateral triangle, as shown by the dashed line in Fig.O 
we call this a 3-point path. We can have n-point paths, 
for any n > 2, which trace out equilateral polygons of n 
vertices. Such paths are interesting because they allow 
us to carry out explicitly the sums in (fl~8"|) and (|2"0|) ; for 
an n-point path, each of the a vectors is equal to one of 
the first n vectors. 

n-point paths for n > 3 are different from 2-point paths 
in that they avoid the center; 3-point paths remain at 
a distance of at least i?i/2 from the central axis, and 
n-point paths for n > 3 remain even farther from the 
center. Another difference is that paths with higher n 
bounce more frequently for the same 9, since they travel 
a shorter distance between bounces. 

In order to test the effects of these differences on our 
rate calculations, we compute 3-point decay and after- 
glow rates by summing (fl~8|) and ([20)) for 3-point paths. 
Our assumption that the particle begins at the center 
of the entrance window gives a valid approximation to 
the decay and afterglow rates only if Tdec^ (2-point) ss 
rdoc, 7 (3-point) and r aft (2-point) « F aft (3-point) . Our 
assumption will be conservative if rd cc , 7 (2-point) > 
rdoc, 7 (3-point) and r a ft (2-point) < r aft (3-point). 

We find that neither the avoidance of the center nor the 
shorter distance between bounces affects the afterglow 
rate very much. The fact that n-point paths for n > 2 
avoid the center is accounted for in the normalization 
factor /geom for these paths. Also, r a f t is dominated by 
paths with small 9, which are more likely to reach the 
PMT. Such paths bounce only a small number of times, 
and the computed r a f t does not change appreciably. The 
3-point afterglow rate is r a f t (3-point) = 7.005 x 10~ 9 Hz 
for the parameters used in Table [TTJ This is within a few 
percent of r aft (2-point) = 6.643 x 10" 9 Hz. 

We also find that the decay rate is somewhat lower 
for 3-point paths, assuming low chameleon masses. For 
the GammeV geometry and the parameters / rc f = 0.53, 
B = 5 Tesla, m off (chamber) = 10~ 4 eV, and /3 7 = 10 12 , 
we find Fdcc.7 (3-point) = 3.800 x 10~ 5 Hz, compared 
to r d cc, 7 (2-point) = 4.905 x 10~ 5 Hz. This is due to 
the fact that 3-point paths bounce more frequently, sup- 
pressing the coherent buildup of photon amplitude over 
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multiple bounces. In our constraints, we use the 2-point 
computation of the decay rate. Since an overestimate of 
rdcc,7 leads to an underestimate of the expected signal 
N(f,(t)T a ft, this approximation is slightly conservative. 

At higher masses, the 3-point paths make a greater 
contribution to the decay rate. As m e g is increased be- 
yond 'ink/L w 10~ 3 eV, destructive interference sup- 
presses oscillation in small-0 paths. Large-0 paths, which 
travel shorter distances between bounces, become more 
important to the computation of Tdoc^- The shortest 
possible distance between bounces for a 2-point path 
is 2i?i, which occurs when 9 = n/2. For m g ~ 
y / Ank/(2Ri) w 10~ 2 eV, oscillation is suppressed by 
destructive interference even for these paths. For a 
3-point path, however, the shortest distance between 
bounces is only Thus rd ec ,7 (3-point) will ex- 

ceed rdec,7 (2-point) at some mass m c g < 10~ 2 eV. Be- 
yond this mass, our approximation of the decay rate by 
Tdoc, 7 (2-point) is no longer conservative. 

GammeV, which probed masses up to 10~ 3 eV, re- 
mained well within the realm of validity of this approx- 
imation; rdoc,7 < Fdcc,7 (2-point) for m c g < 10 -3 eV. 
Thus Fdec,7 (2-point Vwas used to approximate Tdoc^ m 
the analysis of ref. [l7[ in order to provide conservative 
bounds on the chameleon parameter space. We will re- 
visit the issue in Sec. IIVB1 where we shall see that this 
approximation breaks down at the largest masses probed 
by GammeV-CHASE. 



B. GammeV-CHASE 

As discussed in [Tt], HH, the original GammeV 
chameleon search was constrained by four different tech- 
nical limitations: 

1. destructive interference in the L = 6 m mag- 
netic field length suppressed the production of 
chameleons with masses greater than y/ Attcu / L w 
10" 3 eV; 

2. systematic uncertainties in the PMT dominated the 
total error, weakening constraints at low /3 7 ; 

3. the transition between filling the cavity and collect- 
ing afterglow data required to = 1006 sec after the 
laser was turned off, diminishing sensitivity to high 
/3 7 (rapidly decaying chameleons); 

4. the roughing pump in the vacuum system ex- 
hausted to the room, which meant that the lowest- 
density "wall" of the chamber was the P — 1.9 x 
10~ 3 torr intake of the roughing pump. 

GammeV-CHASE improves considerably upon Gam- 
meV by addressing each one of these limitations. 

First, glass windows will divide the magnetic field re- 
gion into partitions of different lengths, 0.3 m, 1.0 m, 
and 4.7 m. These three partitions remove regions of in- 
sensitivity since the chameleon-photon oscillation lengths 



for each partition are not commensurate with each other 
for many multiples. Also, since the partitions are some- 
what shorter than the original 6m cavity, they provide 
some sensitivity to larger mass chameleons, up to a few 
meV. This improvement is especially significant because 
chameleon masses at the dark energy scale 2.4 x 10 -3 eV, 
were inaccessible to the previous experiment. 

Improvements to the optical system will allow 
GammeV-CHASE to push to both lower and higher /3 7 . 
Modulating the PMT signal using a mechanical shutter 
allows the detector noise to be monitored. Since system- 
atic uncertainty in the PMT dark rate was the dominant 
source of noise in GammeV, the sensitivity to low after- 
glow rates will improve by roughly an order of magnitude 
in GammeV-CHASE. Meanwhile, sensitivity to high /3~ 
will improve by three orders of magnitude with the incor- 
poration of two changes. First, a more rapid transition 
will be made between filling the cavity and collecting 
data, shortening the dead time by a factor of more than 
100 and improving sensitivity to larger /3 7 by almost two 
orders of magnitude. Second, data will be collected at 
lower magnetic fields of 1 and 0.2 Tesla, in addition to the 
original field of 5 Tesla. The low magnetic field provides 
an additional order of magnitude in sensitivity because it 
slows conversion of chameleon particles to photons; this 
maintains a detectable population of chameleons while 
the transition to data acquisition occurs. 

Finally, improvements to the pumping system will take 
two forms. First, the vacuum pressure will be reduced by 
approximately three orders of magnitude from ~ 10 -7 
torr to ~ 10~ 10 torr. Second, the vacuum system will 
not exhaust to the room as it did in the original experi- 
ment. These improvements come through the use of ion 
pumps placed at strategic locations on the apparatus, as 
well as cryogenic pumping within the magnetic field re- 
gion; residual gases will freeze to the bore of the magnet. 
The fact that this system does not exhaust to the room 
will mean that chameleons need only bounce from the 
chamber walls, with densities p ~ 1 g/cm 3 , rather than 
the much stronger condition that they bounce before the 
intake of the roughing pump, with p ~ 10 -9 g/cm 3 . Im- 
provements to the pumping system will allow GammeV- 
CHASE to probe chameleons whose masses scale much 
more slowly with density. For chameleons with m e g oc //', 
GammeV-CHASE will probe T) > 0.2. Moreover, the 
improvement in constraints will be qualitative as well 
as quantitative. The large range of masses probed will 
mean that constraints on the photon coupling will be 
only weakly dependent on the matter coupling, and vice 
versa, as we shall see. 

To summarize, there are essentially four important di- 
rections in which the sensitivity of GammeV-CHASE im- 
proves upon GammeV: smaller /3 7 , larger /3 7 , smaller rj, 
and larger m c g (chamber). Figure [5] is a "radar chart" of 
the improvements in these parameters between the two 
experiments. 

We showed in Sec. IIV Al that our choice of initial con- 
ditions, chameleon particles which begin at the center 
of the entrance window, resulted in approximations to 
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FIG. 6: Improvements in sensitivity to /3 7 , r), and m e g with 
GammeV-CHASE. Low f3- 1 improves by roughly an order of 
magnitude. High /3 7 improves by more than three orders of 
magnitude. The parameter r\ improves from 0.8 to nearly 0.2. 
The mass sensitivity improves by a factor of nearly four. 



rdcc,7, r a f t , and the total afterglow flux that were accu- 
rate, and slightly conservative, for GammeV. However, 
at larger chameleon masses, we expect these approxima- 
tions to break down as trajectories with smaller distances 
between bounces contribute more to the decay rate. This 
is because destructive interference suppresses chameleon- 
photon oscillation on segments of the particle trajectory 
that are longer than A-Kk/rri^ s . For example, in GammeV 
and GammeV-CHASE, this length corresponds to 23 cm 
when m c ff = 5 x 10~ 3 eV. 

The n-point trajectories discussed earlier, for n > 3, 
travel smaller distances between bounces. For example, 
3-point trajectories have a minimum distance between 
bounces of Riy/3, compared to 2R\ for the 2-point tra- 
jectories resulting from our initial conditions. At greater 
TO ff, n-point trajectories with greater n will become im- 
portant. By comparing the 2-point and 3-point calcula- 
tions of the decay rate, we estimate the mass at which the 
2-point decay rate is no longer a conservative approxima- 
tion to the total decay rate. Beyond this mass, at which, 
rdoc,7 (3-point) = Tdcc, 7 (2-point), our computation of the 
decay rate becomes increasingly inaccurate. 

Figure [7] shows 2-point and 3-point computations of 
the decay and afterglow rates, for an afterglow exper- 
iment with a single 6 meter partition in the magnetic 
field region. At low mass, it is clear from the figure 
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FIG. 7: rdcc, 7 (thick lines) and T a f t (thin lines) vs. m e ff, for 
2-point and 3-point paths with two different values of £ re f. 
We assume B — 5 Tesla and /3 7 = 10 12 . For simplicity, we 
use a geometry similar to that of GammeV-CHASE but with 
only one L — 6 m partition in the magnetic field region. 



that our approximations rd cc , 7 « F(j 0C 7 (2-point) and 
Taft ~ r a ft (2-point) are excellent. Where the 2-point and 
3-point computations differ, our approximations are con- 
servative; we slightly underestimate the afterglow rate 
and overestimate the decay rate. At larger masses, 
the 2-point decay rates begin to drop, and the 3-point 
rates catch up. Fdec,7 (3-point) exceeds Tdoc^ (2-point) at 
TO c ff ps 4 x 10~ 3 eV for f re f = 0, and at m e ff » 3 x 10~ 3 eV 
for £ re f = 7r; we find similar masses for other values of £ rc f . 

At larger masses, we must compute the decay and 
afterglow rates by averaging ([24]) and (j25|) over all ini- 
tial positions (x, y) on the entrance window. The sym- 
metries that allowed us to sum explicitly the series in 
(TT5)) and (f2"0"|) will no longer be present, making the in- 
tegration more computationally intensive. In this work, 
for the purposes of forecasting GammeV-CHASE con- 
straints, we simply cut off our constraints at the mass 
where Fd 0C ,7 (3-point) first exceeds Fd 0C ,7 (2-point). 



V. CHAMELEON MODELS: CONSTRAINTS 
AND FORECASTS 

A. Model-independent constraints 

The decay rate (|24[) and the afterglow rate ([25|l depend 
on three properties of the chameleon particle: its photon 
coupling /3 7 , its phase shift £ ro f at each wall reflection, 
and its effective mass m c g (chamber) inside the vacuum 
chamber. Constraints on these parameters are model- 
independent in the sense that they do not depend on a 
knowledge of the chameleon potential V(4>). 

Also important is the requirement that chameleon par- 
ticles be contained inside the vacuum chamber, which im- 
plies that the effective mass m e ff (wall) inside the chamber 
walls is greater than the chameleon energy uj — 2.33 eV. 
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FIG. 8: Model-independent 3u constraints from GammeV. 
The shaded region shows the phase-independent constraints 
from [13] , which computed the afterglow rate using only non- 
bouncing trajectories. The solid (red), long dashed (green), 
and short dashed (blue) lines correspond to £ rc f = 0, 7r/3, 
and 7r, respectively. Note that these constraints apply only to 
chameleons which can be contained in the vacuum chamber, 
m e ff (rough) > k, as discussed in Sec. IIVI 



If we assume that the scaling of effective mass with den- 
sity can be approximated by a power law m c g(p) oc p rt 
over the densities of interest, then the smallest value of 
r\ that can be constrained is 
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^?mir 
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(32) 

For j] > ?7 m i n , a larger range of chameleon masses can be 
excluded. 



1. GammeV 

GammeV [ljj constrained the afterglow signal ([3Tj) av- 
eraged over the observation window, which began at a 
time to = 1006 sec after the laser was turned off, and 
had a duration At — 3616 sec. A systematic uncertainty 
of cr r = 12.0 Hz in the PMT dark rate r dark = 115 Hz 
meant that GammeV could rule out to 3cr any chameleon 
model with an average afterglow greater than 36 Hz in 
the time window. 

Figure [8] shows our model-independent constraints 
for GammeV. The shaded region shows GammeV con- 
straints from [171 ] . These are not only model-independent 
but also phase-independent; we considered only non- 
bouncing trajectories when evaluating r a f t , in order to 
remove the effects of phase shifts, and assumed £ rc f = 
when calculating rd cc ,7 in order to overestimate the de- 
cay rate. Since a larger decay rate suppresses the ex- 
pected signal F a f t found in (|3Tj) . this set of assumptions 
is conservative, and applies to all £ re f. 

Meanwhile, the solid and dashed lines in Fig. [5] show 
the regions of parameter space excluded when we con- 
sider bouncing trajectories in the calculation of r a f t and 



assume a specific value for £ re f. Since the decay rate 
rdcc.7 is dominated by paths which bounce hundreds 
times, while paths contributing to the afterglow rate 
r a ft bounce only ~ 1 times, rd 0C ,7 depends much more 
strongly on £ re f than does r a f t . Thus, at high /3 7 , where 
our constraints are limited by rapid chameleon decays, 
the suppression of rd cc ,7 due to nonzero £ rc f extends our 
constraints to stronger photon couplings. On the other 
hand, at low /3 7 , the decay time l/rd ec ,7 is much longer 
than the duration of the experiment, making decays ir- 
relevant. In this regime, the slight suppression of r a f t at 
nonzero £ re f means that the £ re f = constraints are the 
strongest. 

The constrained regions contain several "islands" at 
m j > 10~ 3 eV. These are caused by the zeros of 
the chameleon production probability, given by © with 
t = L. Photons of energy k = 2.33 eV passing through 
a magnetic field region of length L = 6 m cannot pro- 
duce chameleons with masses m e g — \J Airjk/ L, for any 
positive integer j, due to total destructive interference in 
photon-chameleon oscillation. 

The range of rj values in GammeV was severely limited 
by the pumping system used to maintain the vacuum in- 
side the chamber. The lowest-density "wall" was the in- 
take of the roughing pump, with p wa n = 3. Ox 10~ 9 g/cm 3 ; 
chameleons too light to reflect from this intake would 
be pumped out of the chamber. For a reflection phase 
of £ re f = 0, the largest mass probed by GammeV was 
m^ a ^ (chamber) = 1.27 x 10~ 3 eV, giving ry m i„ = 0.76. 



2. GammeV-CHASE 

GammeV-CHASE will monitor the PMT dark rate in 
real time by modulating the afterglow signal from the 
chamber. Furthermore, rapid switching of the PMT will 
allow the experiment to begin data collection at to w 
1 sec, so that GammeV-CHASE will probe chameleon 
theories with much higher decay rates than those ex- 
cluded by GammeV. Here, we estimate the constraints 
that will be obtained by GammeV-CHASE by averaging 
the signal over a time window of At = 20 sec, assum- 
ing that the PMT is modulated with a duty cycle of 0.5. 
The uncertainty in the dark rate Tdaxk ~ 100 Hz over 
this interval will be or ~ 3.16 Hz. 

Forecast constraints from GammeV-CHASE are 
shown in Figure [SI for £ rc f = 0, tt/3, and tt. The shaded 
regions show current constraints from GammeV [l7T | and 
particle accelerators [l||. GammeV-CHASE is expected 
to bridge the gap between GammeV constraints and ac- 
celerator constraints for masses around the dark energy 
scale, TO e ff ~ 10~ 3 eV. 

The constraints shown are limited at low /3 7 by the sen- 
sitivity to low, nearly constant chameleon fluxes. Con- 
sider, for example, £ rc f = 0. From Fig. [3] and the propor- 
tionality Fdcc,7, Taft oc /? 7 , we see that Fd cc ,7 ~ 10~ 6 Hz 
and r a f t ~ 10 -10 Hz when (5 1 — 10 11 and m e g is small. 
The corresponding decay time is much larger 
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FIG. 9: Forecast model-independent 3u constraints from 
GammeV-CHASE. The shaded regions show current con- 
straints from GammeV [TtJ as well as the accelerator con- 
straints of [19). The solid (red), long dashed (green), and 
short dashed (blue) lines correspond to £ rc f = 0, it/3, and it, 
respectively. 



than any time scale associated with the experiment, 
meaning that the number of chameleons in the cham- 
ber, given by (|30|) . remains nearly constant at ivi 8 ^ ss 

EyV^r px ~ 10 n (/Vl0 n ) 2 for the duration of the ex- 
periment. At low /3 7 and m c ff, the total afterglow flux 

reaching the detector is ivj, aft) r aft ~ 10(^ 7 /10 n ) 4 Hz, 
which drops below our sensitivity 3or ~ 9.5 Hz around 
/3 7 = 10 11 . Meanwhile, at high (3 1 , the constraints are 
limited by rapid decays, TJ^ -C to = 1 sec, which re- 
duce the chameleon population before our detector can 
be switched on. Again assuming £ rc f = and low m e ff, we 
see from Fig. Hthat r doCi7 ~ 10- 4 (£/5T) 2 (/3 7 /10 12 ) 2 Hz - 
At the lowest magnetic field used, B — 0.2 Tesla, the 
corresponding decay time rj ec drops below to for /3 7 ~ 
10 , the upper limit of our constrained region. Finally, 
as with GammeV, we see that the £ ro f = constraints 
are strongest at low /3 7 , while constraints at larger phase 
shifts extend to higher /? 7 . 

At high m ff, our constraints are cut off by uncer- 
tainties in the calculation of the decay rate. Our ap- 
proximation that Fdcc.7 < rdcc,7(2-pt) breaks down 
around m c ff « 4 x 10~ 3 eV, where Fd G c,7(3-pt) exceeds 
rdcc,7(2-pt). Had we instead approximated Fdec,7 ~ 
max(Fdcc.7(2-pt), Fdoc,7(3-pt)), our constraints would 
have extended to 6 x 10~ 3 eV. 

Compared with the GammeV constraint plot in Fig.0 
the GammeV-CHASE forecast constraints in Fig.[5]havc 
few sharp features such as islands. This is due to the use 
of partitions in the magnetic field region, as well as data 
runs at multiple magnetic field values. The partitions en- 
sure that there are no zeros in the total chameleon pro- 
duction probability © in the range of masses probed. A 
chameleon whose mass prevents it from being produced 
in the 4.7 m partition, for example, may still be pro- 
duced in the 1.0 m partition or the 30 cm partition. 
Multiple magnetic field values ensure overlap between 
constraints from the different runs. For example, the 



greatest /3 7 probed by the B = 5 Tesla run will be larger 
than the smallest /3 7 probed by the B = 1 Tesla run, 
ensuring a continuous constrained region. Overlaps be- 
tween multiple magnetic fields and multiple partitions al- 
low GammeV-CHASE to smooth over the features that 
are present in the decay rate, and especially in the after- 
glow rate, as in Fig. [3] 

Modifications to the vacuum system, discussed earlier 
in Sec. lIVBl will drastically improve the range of r\ values 
probed. The lowest-density wall will now consist of the 
glass windows, with p wa n ~ 1 g/cm 3 . Inside the magnetic 
field region, the operating temperature of 4 K and pres- 
sure of 10 -10 torr will reduce the matter density inside 
the chamber to p vac = 8 x 10~ 16 g/cm 3 . For chameleons 
with /3 7 <§; /3 m , we can neglect the energy density associ- 
ated with the magnetic field, and approximate m e s oc p^. 
In that limit, ?y m i n = 0.18. For /3 7 ~ f3 m , the magnetic 
field becomes important, since p 1 = 1.11 x 10~ 13 g/cm 3 
is greater than the matter density. In this case, a mass 
scaling m c g oc (p m + p^f 1 is more appropriate, and we 
find 77 min = 0.21. 



B. Power law chameleons 

Consider a chameleon theory with potential V(4>) = 
gcj^ ', for g > and Af > 2. Since /3 7 ^/M P1 , /3 m 0/M P i < 
1 , in order for the theory to be consistent with fifth force 
constraints, the effective potential is given by 



V e s(4>) = g4>- 



(PmPm + /3-yP-y). (33) 



In a bulk medium of constant p m and p 7 , the chameleon 
mass is evaluated at the minimum of this potential, 



(34) 



AT-2 
~ 2 



AfgM Pl 

Note that the mass scales with density as m e ff oc (/3 m pm + 
/3 7 p 7 ) ,? with r] — F° r ^ 4 t neor y we have r\ = 1/3. 

Although r\ grows with Af, it asymptotically approaches 
1/2 at large Af. Thus, the original GammeV experi- 
ment could not exclude any of these power law models. 
GammeV-CHASE, on the other hand, will impose strong 
constraints on power law chameleons, as we will show for 
the model. 

Ref. [HJ finds that the phase shift for power law 
chameleon models, for any real Af, is given by 



tarf - 2 



3Af-2 



Af- 2 



- 1 



(36) 



In particular, for c/> 4 theory we find £ rc f = 2ir, equiva- 
lent to zero. The inverse power law potential V oc 1/0 
has a phase of £ re f = 7r/3. In the limit of large \Af\, 
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Crcf - * 7T. We also find that the exponential potential 
V <x exp(0/M) has £ re f = n. 

Finally, for a power law chameleon with M > 3 
a real number, we must consider the possibility of 
chameleon fragmentation, in which two chameleon parti- 
cles can interact to produce more than two chameleon 
particles. The PMT is sensitive to photons with en- 
ergy uj ss 2.33 eV, the energy of our laser. Repeated 
chameleon fragmentation results in a population of low 
energy chameleons, whose afterglow photons are too low 
in energy to be observed by the PMT. Assuming a frag- 
mentation cross section crf ra g, the fragmentation rate is 
Tfrag = -W^CTfrag/V, where V = 0.030 m 3 is the vol- 
ume of the chamber. Here, represents the number 
of chameleons energetic enough that their afterglow pho- 
tons will be observable by the PMT; although fragmen- 
tation increases the total number of chameleon particles, 
it decreases N^. 

If fragmentation occurs, then the number of de- 
tectable chameleons is given in the production phase 
— T pr < t < of the experiment by 



dt 



= F 7 V pr - r dcc , 7 iV - <7 frag A 2 /F 

(a 2 -6 2 )sinh(^4±^ 



acosh ra^p^ j + 6ginh ( «W*+r pr ) ^ 
where we have defined 

a - 



(37) 
(38) 




b = r doCi7 V/(2o- frag ). 



(39) 
(40) 



In the afterglow phase t > 0, the number N$(t) of 
chameleons and the time average of the afterglow sig- 
nal F a { t (t) = T a ft N^t) over the interval to < t < to + At 
are given by 



N<f,(t) 



2bN <p (0)e- 



26 + A^(o)(l- e - r iec,^) 



(F aft ) = r aft iv (o)/(r£ g ax) Ai) 

nfmax) 



X log 



fra 



(1 - e-rdec.-yAt) 



r 



(max) 
frag 



(i 



(41) 



(42) 



-i(max) 



where the maximum fragmentation rate is If 

A^(0Kag/^. 

Given values of <?, A/", /3 m , and /3 7 , we can determine 
£ re f from (|3"6"|) , and m c ff as a function of the couplings 
from (|35j) . We can compute the afterglow rate and the 
rate of decay to photons for these values of m e ff, /3 7 , and 
£ re f using ((23) and §ZQ . The resulting signal (g2J| can be 
compared to our detection threshold, 3or = 9.5 Hz. In 
this way, constraints on any power law chameleon may 
be obtained. 

Figure [TO] shows the forecast GammeV-CHASE con- 
straints on the <f) A chameleon, with V(4>) — A</> 4 /4!. The 
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FIG. 10: Forecast constraints on the (f> 4 chameleon. Solid 
(red), long dashed (green), medium dashed (blue), and short 
dashed (violet) correspond to A = 10 _1 , 10 -2 , 10 -4 , and 
10~ 6 , respectively; the shaded region shows the accelerator 
constraints of 
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fragmentation interaction cfxj) — > <\>4><\>4> has a cross sec- 
tion <7f rag = afra g A 4 /w 2 , where af rag ~ 1 is a numerical 
factor that cannot be written down in closed form for 
the four-body phase space of outgoing particles; we have 
approximated af rag = 1. We find that fragmentation pre- 
vents us from constraining the model with A = 1 at any 
value of /3 7 , and the models with A = 0.1,0.01 at low 
values of /3 7 , while our constraints for smaller A are not 
limited by fragmentation. 

The plots for A = 10~ 2 , 10~ 4 , and 10~ 6 , shown as 
dashed lines in Fig. [TOJ have a few distinctive features. 
The upper and lower bounds in /? 7 are due to rapid 
chameleon decays and our detector sensitivity, respec- 
tively, as discussed earlier. At low m , we are limited 
by the requirement that chameleon particles be con- 
tained in our vacuum chamber, m e ff(wall) > 2.33 eV. 
Eq. (|35|) with /9 W aii = lg/cm 3 implies that m e ff = 1.6 x 
10" 3 A 1 / 6 /3i /3 eV. For example, for A = 10~ 2 , contain- 
ment requires that j3 m > 3.4 x 10 10 . Meanwhile, in the 
lower right corners of these three plots, corresponding to 
large /? m and small /3 7 , the mass is essentially indepen- 
dent of /3 7 because /3 m p m ^> /3 7 p 7 . Thus the constraints 
in this region resemble the constraints in the lower right 
corner of Fig. [5] where at larger masses we need larger /3 7 
in order to make up for the decline in the afterglow rate. 
When both (3 m and /3 7 are large, the chameleon mass also 
becomes large, pushing us into the regime where 3-point 
decay rates become important and our decay rate com- 
putation becomes unreliable. Thus we have imposed a 
cutoff in the constraints at these masses, as discussed in 
Scc. rvX2l 

The gaps and islands in the constraint plot for A = 0.1 
demand further explanation. We noted earlier that over- 
lapping the constraints from multiple magnetic fields and 
multiple partitions in the magnetic field region helped 
to smooth out the sharp features seen in the afterglow 
and decay rates, Fig. [3] For example, a sudden dip in 
the afterglow due to one partition, if it occurred at a 
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mass where the afterglow due to another partition was 
smooth and sufficiently large, would not show up as a 
feature in our constraint plot. However, rapid fragmen- 
tation at A = 0.1 reduces the afterglow signal from all 
three partitions, at all three B values, so that they no 
longer overlap. We obtain no constraints on the A = 0.1 
models with the B — 5 Tesla run, since that run probes 
low /3 7 , where the afterglow rate is too small to produce 
an observable signal. The gap in the constrained region 
at /3 7 = 1.6 x 10 14 is due to poor overlap between the 
B = 1 Tesla and B = 0.2 Tesla runs. The islands around 
/3 m = 10 16 and /3 7 = 10 15 are due to bumps in r a f t as a 
function of m e ff in the B = 0.2 Tesla run, in a region of 
/3 7 that is well beyond the reach of the B = 1 Tesla run 
due to rapid fragmentation. 

We note that Gamine V-CHASE constraints on the </> 4 
chameleon will be at very different couplings from the 
constraints of laboratory fifth force searches [HI, |36| . 
Furthermore, a 4 chameleon massive enough to be 
contained in the vacuum chamber, m e g(p wa i\) > uj = 
2.33 eV, will have a mass in the Galaxy of m e g (Galaxy) > 
10~ 8 eV. Since this is much larger than the plasma 
frequency cjp « 10~ n eV of gas in the galaxy, the 
chameleon-photon mixing angle will be strongly sup- 
pressed. Thus, GammeV-CHASE will also probe </> 4 
chameleon models different from those ruled out by as- 
trophysical dimming constraints. 



C. Chameleon dark energy 

1. Inverse power law potentials 

First, we consider an inverse power law potential with 
a constant term, 



V{4>) = Mi 



1/4 



(43) 



where n > 0, k > 0, and M A = = 2.4 x 10~ 3 eV. The 
model with k = 1 is the most economical, in the sense 
that the constant and power law terms both see the same 
mass scale, Ma; however, we allow for the possibility that 
k ^ 1. Eq. (J3SJ) for J\f = —n implies that the phase shift 
at a wall reflection is 



£ref(™) 



rm 



(44) 



The bulk field value m i n , where dV c g {4> m in) / d(j> = 0, and 
the effective mass r?i 2 ff = d 2 V e s (0 m in)/<9</> 2 are given by 



0n „„ = Mh ( ;" M 'f* ) 

\PmPm + P 7 /0 7 / 

m eS = m (n + l)M A l ^ 



(45) 



(46) 




The mass scaling with density is given by r/ — (n + 
2)/(2n + 2), which is 3/4 for n = 1 and 2/3 for n = 2. 



le+11 



10000 1e+06 1e+08 1e+10 1e+12 1e+14 1e+16 
matter coupling pi m 

FIG. 11: Forecast Gamme V-CHASE constraints on 
chameleon dark energy with an inverse power law potential 
(|43[) . Thick and thin lines refer to k = 1 and 0.1, respectively. 
Solid (red) and dashed (green) lines refer to n = 1 and n = 2, 
respectively. 



Only models with n somewhat less than 1 were accessible 
to GammeV, whereas all n will be probed by Gamme V- 
CHASE. 

At the cosmological matter density p cos « 2.5 x 
10~ 30 g/cm 3 , the background field value for the n = 1 
model is <j> min = 1.7 x 10~ 15 (k/ 'Pn^^Mpi, and the 
chameleon Compton wavelength is ~ 100 pc x 

(K//3 m )( 2n+2 )/(" +2 ). The Compton wavelength is tiny on 
cosmological scales, so the field remains close to the min- 
imum of its potential. Most of the chameleon energy 
density will come from the constant term M\ in the po- 
tential. For example, for n = 1, (3 m = 10 15 , and n = 1, 
the field-dependent term k(M\/ '</> m i n ) n ~ 10~ 8 <§; 1, and 
this term is even smaller for larger n and lower /3 m . The 
field dependent term does not exceed the constant term 
until the matter density becomes of order 10 -16 g/cm 3 , 
or about fourteen orders of magnitude greater than the 
current cosmological background density. For all cosmo- 
logical purposes, such a chameleon model behaves like a 
cosmological constant. Evolution in its equation of state 
will be undetectable. 

Since <f> approaches zero at large densities, the total 
fractional change in the fine structure constant, between 
large densities and the cosmological background density 
p cos , is given by Aoem/ ^em = /3 7 ^ m in(pco S )/M P1 . This 

is Aa EM /a E M = 2.4 x 10- 15 /3 7 (K//3 m ) 1/2 for n = 1, 
and even smaller for larger n. Inverse power law mod- 
els where /3 7 is sufficiently greater than /3 m can be con- 
strained by cosmological bounds on the variation in 
qem [13, However, we note that chameleon mod- 

els predict a density-dependence of cvem, rather than 
a time-dependence. Since astrophysical measurements 
of oem are made in overdense regions of the universe, 
Pm > Pcos, the actual variation in q E m expected from 
chameleon theories will be smaller than the value given 
above. Similarly, laboratory searches [39( for chem vari- 
ations carried out under conditions of constant density 
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will not be sensitive to chameleons, though it may be 
possible to derive a constraint by comparing several «em 
measurements such as conducted at slightly different 
densities. Laboratory and astrophysical constraints on 
chameleons from variations in oem will require further 
analysis that is beyond the scope of this paper. 

Photons from extragalactic sources, passing through 
the galactic magnetic field, could potentially oscillate 
into chameleon particles, causing the sources to appear 
dimmed, and possibly polarized [2(| [H|. Constraints 
on astrophysical dimming have ruled out a wide range 
of photon couplings. However, these constraints only 
apply to chameleons with masses much less than the 
Galactic plasma frequency ujp « 10~ n eV; larger masses 
strongly suppress chameleon-photon mixing in the low 
(B ~ 10 -10 Tesla) magnetic fields found in the Galaxy. 
Assuming a Galactic matter density of ~ 10 7 p cos , the 
chameleon mass becomes greater than wp for (3 m > 10 12 
for n = 1 and k = 1. Larger values of n and smaller 
values of k cause m e s to exceed wp for even lower values 
of f3 m ; for n = 2 and k = 0.1, we find that (3 m = 10 10 
is sufficient. Furthermore, astrophysical dimming con- 
straints can be evaded altogether if the chameleon has 
a bare mass greater than wp, in addition to the mass it 
acquires from its matter coupling. 

GammeV-CHASE will complement current con- 
straints on inverse power law chameleons. Figure 1111 
shows forecast GammeV-CHASE constraints for n = 1, 
2, and k = 1, 0.1. The constraints are limited at low 
/3 7 by low afterglow rates, and at high /? 7 by rapid 
chameleon decays. At low matter couplings, (3 m < 10 4 , 
the chameleon mass in the walls of the vacuum cham- 
ber drops below their energy u>, preventing the cham- 
ber from trapping the chameleon particles. At large 
matter couplings, f3 m > 10 15 , the constraints are lim- 
ited by destructive interference due to large chameleon 
masses, and by the breakdown of our approximation 

rdcc, 7 (2-point) w r d0Ci7 . 



2. Exponential potentials 

Next, we study chameleon dark energy with an expo- 
nential potential, 



v(4>) = Mi 



K(b 



(47) 



where, once again, Ma = 2.4 x 10~ 3 eV, and k > 0. 
The constant term M^ ensures that <j> behaves as a dark 
energy; without it, the energy density of the field would 
scale with the background matter density, and acceler- 
ated expansion would not take place. For this potential, 
we find, using the techniques of [18J, that the phase shift 
at a wall reflection is £ rc f = 7T. 

The bulk field value and the corresponding effective 
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FIG. 12: Forecast GammeV-CHASE constraints on 
chameleon dark energy with an exponential potential. Solid 
(red), long dashed (green), medium dashed (blue), and short 
dashed (violet) lines refer to k = 1, 10" 1 , 10 -2 , and 10~ 3 , 
respectively. 



V/, 



mass are given by 



Ma 



log 



kM?M P i 



'off 



W (AnPm + PjPy) 

M A Mpi 



(48) 
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The mass scales as the square root of the density, making 
exponential potentials inaccessible to GammeV, but well 
within the abilities of GammeV-CHASE. 

The field value changes very little over a large range 
of densities. The total variation in the log term in m ; n , 
over a range of densities from p cos to laboratory densities 
~ lg/cm 3 , is less than 100. Thus the total change in </> m i n 
between these densities is of the order of 10 _28 k~ 1 Mpi. 
Such a small change in m i n means that cosmological 
variations in a^M will not be observable for /3 7 < 10 23 , 
assuming that k is of order unity [53, HH • 

At the cosmological density, the chameleon Compton 
wavelength is = 4.7 x 10~ 6 pc x (K/3 m ) -1 / 2 , so the 
field remains close to <^ m in(Pcos) on cosmological scales. 
As with inverse power law chameleons in the previous 
section, the field-dependent term in the potential is small 
compared to the constant term, exp(— K(j)/M\) = 3.2 x 
10~ 31 (3 m /n <C 1, so the field behaves as a cosmological 
constant. 

For the exponential chameleon, constraints from astro- 
physical dimming do not overlap with constraints from 
GammeV-CHASE. Any chameleon massive enough to 
be trapped in the GammeV chamber will have a mass 
greater than 10 -11 eV at galactic densities, and dimming 
will be suppressed. 

Figure [T^] shows forecast GammeV-CHASE con- 
straints on chameleon dark energy with an exponential 
potential. As in the inverse power law case, our con- 
straints are limited at low /3 7 by a low afterglow rate, 
at high (3ry by rapid chameleon decays, at low /3 m by 
chameleons too light to be trapped in our chamber, 
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and at high /3 7 by destructive interference in chameleon- 
photon oscillation and the breakdown of our approxima- 
tion that rdcc,7(2-point) rs r<j eCj7 . 

VI. CONCLUSION 

GammeV and GammeV-CHASE are chameleon after- 
glow experiments which probe the coupling between pho- 
tons and chameleon particles that are trapped inside a 
vacuum chamber. Using the physics of chameleon-photon 
oscillation, we have computed the rate rdec,7 at which a 
chameleon converts to a photon inside the vacuum cham- 
ber of such an experiment, and the rate r a f t of pro- 
duction of photons which escape the chamber to reach 
an external detector. We have computed these rates as 
functions of the chameleon-photon coupling /3 7 , the ef- 
fective chameleon mass m c ff (chamber) inside the vacuum 
chamber, and the phase £ re f that is introduced between 
chameleons and photons as they reflect from the walls of 
the vacuum chamber. The decay and afterglow rates for 
the GammeV geometry are plotted as functions of m e g 
in Fig. [3] and £ rc f in Fig. |U With these rates, we have 
computed the number of chameleons in the chamber and 
the expected afterglow flux in our detector, as a function 
of /3 7 , m c ff (chamber), and £ rc f. 

In a previous result, GammeV [ijj used the decay 
rate computed above, as well as the afterglow rate due 
to non-bouncing trajectories, to exclude a region of the 
chameleon parameter space. Here, we have included 
bouncing trajectories in our afterglow calculation, result- 
ing in constraints that are less conservative and that de- 
pend on the phase shift £ re f. These model- independent 
constraints are shown in Fig. [5J For £ rc f = and small 
masses, the new constraints extend from half of the mini- 
mum photon coupling excluded by [ItJ to twice the max- 
imum coupling excluded by that reference. For nonzero 
phases, additional suppression of the chameleon decay 
rate to photons allows us to push to even higher /3 7 , 
excluding photon couplings ten times as high as those 
ruled out by [l7j . However, since GammeV was only able 
to probe chameleon models whose masses scaled rapidly 
with density, even this new analysis does not allow us to 
constrain the most common types of chameleon poten- 
tials. 

GammeV-CHASE is an improved version of GammeV 
that is expected to take data in the winter of 2009-2010. 
In Fig. [9] we forecast the model-independent chameleon 
constraints that will be achieved by GammeV-CHASE. 
These constraints span many orders of magnitude in /3 7 , 
bridging the gap between the constraints from Gam- 
meV and from particle accelerators. We show that 
the improved control of PMT systematics will allow 
GammeV-CHASE to probe smaller couplings than Gam- 
meV. Rapid PMT switching and multiple runs at lower 
magnetic fields will allow us to constrain /3 7 as high as 
10 16 . Our forecasts show that GammeV-CHASE will be 
able to probe a large range of couplings even at masses as 
high as the dark energy mass scale, 2.4 x 10~ 3 eV. Fur- 



thermore, improvements to the pumping system allow 
GammeV-CHASE to probe many commonly used and 
well understood chameleon potentials, including quartic 
potentials, as well as inverse power law and exponential 
potentials that could explain the observed cosmic accel- 
eration. 

Quartic chameleons have been discussed extensively in 
the literature [31], 

H HE S3 • Torsion pendulum experi- 
ments such as |33[ have ruled out matter couplings up to 
[3 m = 1, but arc insensitive to chameleons with stronger 
matter couplings. It is precisely these chameleons which 
will be trapped in the GammeV-CHASE chamber. We 
have shown that GammeV-CHASE will complement 
laboratory searches for quartic chameleons by probing 
strongly coupled chameleons, as shown in Fig. [TDJ For 
self-couplings a few orders of magnitude smaller than 
unity, constraints on quartic chameleons will span seven 
orders of magnitude in f3 m and five orders of magnitude in 
/3 7 , extending from the upper bound /3 7 ~ 10 16 provided 
by particle accelerators down to (3 1 ~ 10 11 . GammeV- 
CHASE is also complementary to searches for a vary- 
ing fine structure constant. If we extend a simple model 
of ckem variation such as [42| to a chameleon theory, by 
adding a 4 potential and a Yukawa matter coupling with 
coupling constant proportional to mass, then constraints 
from GammeV-CHASE will be more powerful than those 
from laboratory or cosmological tests in a large portion 
of the parameter space. 

Chameleon dark energy, with runaway potentials of 
the form V(cj>) = Af|/(0/M A ) and M A = 2.4 x 10" 3 eV, 
have also been studied in the literature [3]. In such 
potentials, /(</>/M A ) — > 1 as 4> runs off to large values, 
corresponding to low matter densities such as the cur- 
rent cosmological background density. We have applied 
our forecast GammeV-CHASE constraints to chameleon 
dark energy models with inverse power law potentials 
V(4>) — Af£[l + (Ma/0)' 1 ] as well as exponential poten- 
tials V{4>) = M£[l + exp(-0/M A )]. Figure HU shows the 
constraints that GammeV-CHASE will be able to place 
on inverse power law chameleon dark energy. For the 
simplest model, n = 1, we will be able to probe matter 
couplings ranging from /? m < 10 5 to /? m < 10 16 and from 
Pi ^ 10 n all the way up to the accelerator constraints, 
/3 7 ~ 10 16 . At the largest matter ouplings, (3 m > 10 12 , 
we will probe parameters inaccessible to astrophysical 
dimming constraints and bounds on variations in the 
fine structure constant. GammeV-CHASE constraints 
on dark energy with an exponential potential, shown in 
Fig. 1121 will be similarly powerful, covering matter cou- 
plings in the range 10 6 < /3 m < 10 16 . Furthermore, 
this range of matter couplings is completely inaccessi- 
ble to astrophysical dimming constraints, which apply to 
chameleons with low masses in the galaxy, and bounds 
from ftEM variation, applicable to models with large vari- 
ations in 0. Thus we have shown that GammeV-CHASE 
will probe large ranges of previously unexplored param- 
eter space for the simplest models of chameleon dark en- 
ergy. 
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